Optimal operating conditions of an entangling 
two-transmon gate 

Antonio D'Arrigo 

E-mail; adarrigo@dmfci.unict.it 

Dipartimcnto di Fisica c Astronomia, Univcrsita di Catania, c/o Vialc A. Doria 6, 
Ed. 10, 95125 Catania, Italy and CNR - IMM - MATIS c/o DFA Via Santa Sofia 64, 
95123 Catania, Italy 

jL( ■ Elisabetta Paladino 

a; 

E-mail: epaladino@dinfci.unict.it 

Dipartimcnto di Fisica c Astronomia, Univcrsita di Catania, c/o Vialc A. Doria 6, 
Ed. 10, 95125 Catania, Italy and CNR - IMM - MATIS c/o DFA Via Santa Sofia 64, 
95123; Catania, Italy 

Abstract. Wc identify optimal operating conditions of an entangling two-qubit gate 
realized by a capacitive coupling of two superconducting charge qubits in a transmission 
line resonator (the so called " transmons" ) . We demonstrate that the sensitivity of the 
optimized gate to 1// flux and critical current noise is suppressed to leading order. 
The procedure only requires a preliminary estimate of the 1// noise amplitudes. No 
additional control or bias line beyond those used for the manipulation of individual 
qubits are needed. The proposed optimization is effective also in the presence of 
relaxation processes and of spontaneous emission through the resonator (Purcell effect). 



PACS numbers: 03.67.Lx, 85.25.-j,03.65.Yz 
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1. Introduction 

Superconducting circuits are a promising technology for the realization of quantum 
information on a solid state platform. Several types of qubits [1] have been developed 
realizing high fidelity single qubit operations [21 [3]. Rapid progress has also been 
made towards the realization of robust and scalable universal two-qubit gates [U Ej E]. 
The circuit quantum electrodynamics (cQED) [7] architecture demonstrated to be 
particularly promising for scalable quantum information. In this scheme highly 
entangled two [HI [9] and three qubits [10] have been generated and simple quantum 
algorithms have been demonstrated [TTl [T2] . 

The coherence times of the present generation of devices (~ /is) are about three 
orders of magnitudes larger than the first implementations. A relevant step further 
toward this enhancement has been the elimination of linear sensitivity to low-frequency 
(1//) noise by operating qubits at "optimal" working points. After the first "sweet 
spot" operation demonstrated in Ref. [3], a further boost of qubit performances has 
been achieved in a cQED design named "transmon" [13], which is almost insensitive to 
the detrimental effect of 1// charge noise [H] at the price of reduced anharmonicity. 
However, cQED architectures share with other implementations the presence of 1// 
flux noise whose amplitude has a characteristic order of magnitude [15], and of 1// 
critical current noise [13]. Together with relaxation processes due to quantum noise, 
dephasing due to 1/f flux and critical current noise still limits the time scales over 
which phase coherence and entanglement are preserved. In fact, further improvement 
of the coherence times at least of one order of magnitude would be required to reach 
the level for practical quantum error correction [16]. Recently in a new circuit-QED 
architecture employing a three-dimensional resonator the error correction threshold has 
been approached [17]. "Optimization" is thus a key- word of the present generation of 
superconducting nano-circuits. Clever circuit design and optimal tuning of multi-qubit 
architectures, supplemented by the use of improved materials, are two complementary 
strategies currently exploited to address this problem. 

A major question currently unsolved is establishing the best strategy to maintain 
long-enough a sufficient degree of entanglement. In the present article we address this 
issue considering a universal two-qubit gate realized by a fixed capacitive coupling of two 
transmons in a cQED architecture. The implementation of this scheme has been recently 
reported in Ref. [9] where a a/I — SWAP operation with individual single-shot non- 
destructive readout [18] and gate fidelity of 90%, partly limited by qubit decoherence, 
has been demonstrated. A similar system has been studied theoretically in [19l |20] . 
Here we identify "optimal" [21] operating conditions of a transmons a/I — SWAP gate 
taking into account the multi-level nature of the nano-circuit. We find that an "optimal 
coupling" exists where the leading order effects of 1// flux and critical current noise 
are eliminated. The amount of preserved entanglement is quantifled by the concurrence 
between the two transmons, C(t), which we evaluate in analytic form. The efficiency of 
the "optimal coupling" is demonstrated by the fact that, for typical 1/f noise spectra 
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measured in superconducting nanocircuits, the concurrence is predicted to decay on 
a time scale [22] Tg*^^"^^ > 300 /is (in the absence of other decay mechanisms). In 
addition, C{t) may attain values [23] guaranteeing violation of a Bell inequality until 
~ 80 /is and the gate fidelity is 99% up to ~ 20yus. Finally, we demonstrate that 
the optimization is effective also in the presence of relaxation processes due to flux 
quantum noise. Similarly to other cQED systems [21], the gate efficiency can be limited 
by spontaneous emission through the resonator. This limitation is likely to be overcome 
by suitable Purcell filters or protected designs [25]. The optimization proposed in the 
present article can further improve the considerable performance of cQED two-qubit 
gates based on cavity-mediated interaction [T2l [TO] or on tunable effective interaction 
with microwave control [26]. Remarkably, here effective elimination of omnipresent 1// 
noise sources is achieved even if one qubit does not operate at optimal bias and without 
additional controls or bias lines beyond those used for the manipulation of individual 
qubits, an important feature for scalability. 

2. Universal two-transmon gate 

We consider two transmons with a fixed capacitive coupling, each qubit being embedded 
in its superconducting resonator used for control and bit-wise readout [9l [18]. The 
interaction is effectively switched on/off by dynamically changing the qubits detuning 
using single qubits control lines. For this reason one of the qubits does not operate at 
its sweet spot. In figure [1] (a) we report the circuit diagram of the considered system. 
Each transmon, denoted by the subfix a = 1,2, consists of a Cooper-Pair-Box (CPB) 
characterised by the charging energy Eca and Josephson energy Ej^ = Ej^cos{(j)a), 
tunable via the magnetic flux threading the superconducting loop, (pa = tt^q/^o? 
($0 is the flux quantum). In the circuit-QED scheme each CPB is embedded in a 
transmission line resonator whose relevant mode is modeled as a LC oscillator |13j . 
Thus the Hamiltonian of transmon a consists of the CPB Hamiltonian plus the dipole- 
like interaction with the LC oscillator [13] 

Ha = Eca{qa " qx,af " Eja{(f)a) COS (^q + Ura a^tta + WaGVa Qa (Oq + «!) , (1) 

where phase, (pa and charge, g^, are conjugate variables, [paAc] = ^- The resonator 
energy is Ura = ^/ V^aCa and a]j (a^) creates (annihilates) one photon in the 
transmission line {h = 1). Va = \/ Uroi/'^Ca is the root-mean-square voltage of the 
oscillator and /3q, = Cga/Cj^a is the ratio between the gate capacitance coupling the 
CPB to the local mode and the CPB total capacitance. 

The transmon operates at Eca ^ Ej^. Under these conditions, values of the 
phases close to zero are most favored. This motivates the neglect of the periodic 
boundary condition on the phases and the expansion of the cosine in Eq. ([T]). Within this 
approximation, the offset charge q^^a can be eliminated via a gauge transformation [13] . 
Of course, the perturbative scheme cannot capture the non-vanishing charge dispersion 
of the transmon [27| [28] . In particular, the exponential decrease of the charge dispersion 
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Figure 1. (a) Circuit diagram of the two transnion gate with capacitive coupling 
energy Ecc = {2e)^CT /CsiC^2, where 1/Ct = l/Cc + 1/Csi + 1/Cs2 and = 
Cja + Cga + Cbo- (b) Schematic level structure. 



with ^y Ej^/Eca for Ej^/Eca ^ 1 only results from the exact diagonalization of 
the CPB Hamiltonian in the phase basis [13]. It leads to exponential suppression of 
sensitivity to low-frequency (1//) fluctuactions of the offset charge [2]. Here we rely 
on this well established result and eliminate q^.a from the outset. Expanding the cosine 
in Eq. ([T]) up to fourth order, the CPB Hamiltonian can be cast in the form of a weakly 
anharmonic oscillator (Duffin oscillator) 

= bib^ - {Ecaim{ba + blY , (2) 

where the bosonic operators b^, b]^ are related to the charge operator via qa = 
-i{Ejj2Ecafl\b^ - bi)/V2 and we put = = ^2EcaEjMa)- The 

two lowest eigenenergies of "H^ identify the transmon-a qubit levels. Their splitting 
is Qa = — Eca/4: and it can be tuned by changing the magnetic flux The flux 
"sweet-spot" is at 0q, = [3l [13]. 

The capacitive coupling between the CPBs, Ecc{Qi — (1x,i){Q2 — Qx,2), adds to Tia 
leading to the Hamiltonian 

'H = 'H? + n^ + ^ (&1 - &I)(&2 - bl) , (3) 

where Ecc = Ecc{EjiEj2/AEciEc2Y^'^ is the effective coupling depending on the 
control parameters 0^ via the Josephson energies. Note that fluctuations of the 
magnetic fluxes affect the effective coupling between the qubits. Typical values [9] 
are Eca ~ 1 GHz and E^j^ - 20 - 30 GHz, leading to Ecc = 10"^ - lO'^GHz. 

The coupled transmons eigenenergies and eigenstates are conveniently obtained by 
treating in perturbation theory with respect to J2a b^^ba both the anharmonic terms 
and the capacitive interaction included in 

v = ^{b^- bi)ib, - 4) - E ^{b. + bir . (4) 

a 

The level structure is schematically show in flgure [1] (b). The splitting in the 
subspace where the SWAP operation takes place (in short "SWAP splitting") reads 

aj+„ = ■\J (ill — 1^2)^ + E'^c ^^'^ ^1^^ corresponding eigenstates spanning the "SWAP 
subspace" are |-) = - sin(r//2)|01) + cos(r//2)|10) and |+) = cos(r//2)|01)+sin(r//2)|10). 
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where tan 77 = Ecc/iP-i — ^2) and \a,h) = |a)i|6)2 are eigenstates of '^^D.ab'l^ba, (a, 
b G {0, 1}). The interaction is effectively switched on by tuning the single-qubit energy 
spacing to mutual resonance. The resonance condition is realized by tuning the flux bias 
until Cli = CI2, displacing one qubit from the sweet spot at 0^ = 0. In the following we 
suppose that 0i = and 02 7^ 0. Under resonance conditions the \/i — SWAP operation 
|01) — )■ \^e) = [|01) — z|10)]/a/2 is reahzed by free evolution for a time = 7r/2a;+_ 
starting from a factorized initial state in the " SWAP subspace" . 



3. Optimal operating conditions: reduction of 1// noise effects 

Since the two qubits do not operate at the same working point, the dominant source 
of dephasing is different for the two transmons. In particular, first order fiuctuactions 
of the transmon splittings are due to 1/f critical current noise for transmon 1 and 
to 1/f flux noise for transmon 2 [13]. These fluctuations can be treated in the 
adiabatic and longitudinal approximation [29] by replacing Eja with Eja{l + Xait)). 
Here Xa{t) represent stochastic fluctuations of the dimensionless critical current Xi{t) = 
SIci(t) = AIci(t)/Ici, and of the flux $2, X2{t) = tan(02) ^02(^)- The leading order 
effect of adiabatic noise is defocusing, expressed by the "static path" or static noise 
approximation (SPA) [291 EHj describing the average of signals oscillating at randomly 



distributed effective frequencies (see Appendix A for the validity regimes of the SPA 



in the present problem). It is obtained by replacing Xa{t) with statistically distributed 
values Xa- In the SPA the coherence between the states |±) is 

(p+_(i)) = P+40) e-^-+-* (e-^-+-*) , (5) 

where p{t) is the two-qubit density matrix and (...) indicates the average over the 
fluctuations Xa- Here we assume that they are uncorrelated random variables with 
Gaussian distribution, zero mean and standard deviations proportional to the 
amplitude of the 1/f spectrum, S^J {uj) = 7r'El^[\n{'yMa/'yma) '^]~'^ {ima and 7^0 are the 
low and the high frequency cut-offs of the 1/f region). As demonstrated in Refs. [2X1122] 
the optimal operating condition is obtained imposing a minimum of the variance of the 
stochastic SWAP splitting, T? = {<jj\J) — (cj+_)^. This is simply understood considering 
the short-times expansion |(e^*'^'^+"*)| ~ ^/Y^^JT/Vf , implying defocusing suppression 
when E is minimal. Expanding around the fixed working point we get 

a ^ ' a, 13 ' 

where all derivatives are evaluated at x^ = 0. At resonance we find duo^^ / dxa = Eqc/^i 
d^u+^/dxl = -3Ecc/lQ + ^Vi^Ecc), d'io+^/dx,dx2 = -Ecc/IQ - ^Vi^Ecc), 
where we put = Eca = Ec- The variance Eq. (|6l) is non-monotonic in the 
coupling energy (figure [21 (a)) and its minimum depends on the noise variances S^^. 
For typical values of the amplitudes of 1// flux and critical current noise [32] the 
dominant effect is due to flux noise, S^.^ ^ S^.^ and the optimal coupling is found at 
E^P^ ^ 2Ec{T.xJ V^f'^ ■ Note that, since Ecc depends on Xq, (via Eja), the differential 
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Figure 2. Panel(a): SWAP-splitting variance Eq.® as a function of Ecc/Ec- 
The value of the minimum of S at the optimal point E'^^ = 2Ec{'£'x2/ V^Y^"^ = 
1.68 • 10~^Ec is one order of magnitude smaller than at Ecc = lO~^Ec- Panel (b): 

Dispersion branch Slu^ {x2, xi = 0)/ Ec for \x2\ < SE-c^. The black line is for a generic 

coupling Ecc ~ 10^^i?Ci the red line is for the optimal coupling E'^^. Parameters 
are Eca = IGHz, E^j^ = 30Eca with ^2 = 0.64 and T.^,^ = 10"^. 



dispersion duj+-/dxa at = is non vanishing unless the couphng is switched off. 
The condition of minimal variance effectively identifies an "optimal" dispersion leading 
to minimal defocusing, see figure |2t^b). 

In addition we observe that, since E'^^ depends on Ec but not on the Josephson 
energy, the optimized SWAP frequency, 0;°"^ ~ -^cc' engineered by appropriately 

fixing (within the experimental tolerances) the ratios Ej^/Eca- This recipe can be 
conveniently applied even if an independent estimate of the flux noise amplitude, S^^j , 
for the specific setup is not available. In fact, the variance of the stochastic SWAP 
splitting, S^, depends very smoothly on Ecc (figure |2] (a)) allowing a practical estimate 
of E^Q based on the characteristic value of S^-j observed in different flux and phase 
qubits. Alternatively, if different devices can be fabricated, one should select the sample 
with the ratio Ej^/Eca taking the right value for the given noise level of that particular 
device. 

The effectiveness against defocusing of operating at the optimal coupling is revealed 
by the concurrence [33], which we evaluate in the SPA. We assume the system is prepared 
in the state |01) and freely evolves. In the adiabatic approximation populations are 
constant thus C{t) ~ 2 |Im(p+_(t))|. Evaluating the integral ([5]) we obtain 



A measure of the entanglement preservation is the "SWAP decay time" [22] defined 

by the condition | Cspa (^'2*^^'^^ ) | = e~^. At the optimal coupling is one order 

of magnitude larger than for a generic coupling, assuming remarkable values up to 
j^*swAP _ 4QQ^g g^g^^g ^^^Yi increasing E^Ec, figure [Sj^a). In addition, a 99% fidelity 
to the Bell state \ipent) = [\0l) + \l0)]/\/2 is maintained up to Tjrgg ^ 20 /xs, about 4 times 
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Figure 3. Tf'^^'^ (panel a) and Tjc-gg (panel b) as a function of Ey Ec for 
Eca = IGHz. The variances E^.^ = 5 • 10~^, T^^^ = 10"'* correspond to typical 
values of 1// critical current and flux noise [35] rescaled to the present setup. Black 
lines are obtained for Ecc = ^^~^Ec, red lines correspond to E'^^ = 1.68 • IQ^'^Ec- 

longer than for a generic coupling, figure |3l^b). These results elucidate the capability of 
the proposed operating condition to drastically reduce defocusing due to 1// flux and 
critical current noise. On the other hand, energy relaxation processes are expected to 
limit the gate fidelity and the qubit relaxation times of the considered architecture [9] . 
In the following Section we discuss the robustness of the optimal coupling condition to 
relaxation processes. 

4. Optimal operating conditions: robustness to relaxation processes 

We now discuss the robustness of the above optimization against relaxation processes 
due to flux noise and to spontaneous emission through the resonator. Flux quantum 
noise is due to the external magnetic flux bias through a mutual inductance M [13] and 
it enters the Josephson energies in T/d and in Eqc- It is included by adding to "H, Eq. 
(I3|), the terms 

= 5^f^a&I&aX, - ^ (^1 - - h\) J2^o. • (8) 

a a 

For the transmon at the flux sweet spot it is Xi = 50^/2, for transmon 2 instead 
X2 = tan 02(502, where 6(f)a are quantized phase fluctuations. A'H conserves the parity 
of the total number of the two transmons excitations. Thus it does not connect 
the states |±) to the ground state which, to the flrst order in V takes the form 
\g) oc |00) + aii|ll) + X]ij=i 2 '^ojlO, 2j) + aio\2i,0). Disregarding thermal excitation 
processes to higher energy states, the only effect of flux quantum noise is inside the 
bi-dimensional subspace {|±)}. By solving a Bloch-Redfield master equation [M] 
relaxation and decoherence times in the SWAP subspace are given by the usual relation 
T|WAP = 2TfWAP = {±J2M^V+{-'^rcosr]Ecc/^afS,^{u+^)}-' (apure dephasing 
term oc 5'a;^(0) is disregarded with respect to defocusing due to 1// noise). 
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Figure 4. Envelope of the concurrence in the presence of 1// flux and critical current 
noise (as in figure[3]) and flux quantum noise on qubit 2 with spectrum ('^) — f O^^w. 
Black line is for Ecc = lO^^-Bci; the red line is for E°^^ = l.QS-lQ^'^Eci. The dashed 
gray line marks the value C = 1/V2, the dotted hue C{Tf^^^) = e~^. 



The main contribution comes from linear phase fluctuactions of the transmon 
displaced from the sweet spot, X2. At low temperatures fc^T <C flux quantum noise is 
Sxii.^) ~ (^^2tan(/)2)^(7rM/$o)^2cu/i?. For typical parameters we estimate T^"^^^ ^ 30s 
at optimal couphng (M = lAO^o/A, R ~ 50f2) [13]. Thus the efficiency of the optimized 
gate on the SWAP time scale is not limited by relaxation processes due to flux noise. 
This is illustrated in figure H] where we plot the envelope of the concurrence, C{t), which 
in the presence of 1// noise and fiux quantum noise reads 



C{t) ^ {[sinr/(p++(t) - p__(t)) + 2cosr/Re{p+_(t)}]2 + Almip^^t)}'}'/' . (9) 



Here the population difference in the SWAP subspace is p++(t) — p (t) = 

(cosr/ — 5eq)e~^^'^^^^^ + Seq, with 5eq the thermal equilibrium value, and p+-(t) = 
(p+_(t))e~*/^2 . We observe that for optimal coupling Bell inequality violation, 
guaranteed until [35] C{t) > 2"^/^, occurs for times ~ 75 ps, much longer than for 
generic coupling. 

In the cQED architecture each transmon is dispersively coupled to a resonator 
used for control and readout. An important mechanism for Ti processes is spontaneous 
emission through the resonator (Purcell effect) [131 IM]- If each transmon operates at 
positive resonator-transmon detunings, Aq = Ura — ~ 500MHz ^ (7a ~ 50 MHz, 
where ga is the transmon-resonator coupling strength, the spontaneous emission rate of 
the coupled transmons is due to the "single-mode" Purcell effect [241 IlH]- To evaluate 
it we rewrite H, Eq. ([3]), in the basis of its perturbative eigenstates and perform the 
rotating wave approximation eliminating terms describing the simultaneous excitation 
(de-excitation) of one resonator and the coupled-transmons system. The restriction to 
the subspace (l^f), |cr = ±)} reads 




(10) 
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where ga^^j = i\/2l3aeVa{Eja/2EcaY^^{g\ipa — b\)\a). The eigenstates of (fTUi) are 
obtained by treating the last term in first order perturbation theory. The ground 
state is unmodified and reads |(7,0i,02), where \ma) are Fock states of the a-th 
resonator, G N. The corrections to the states |cr, mi,m2) read | a, mi, 7712) = 
gia\/rn.i+i _| mi + 1, ^2)+ ^2<T\/m2+i | m2 + l), where and E„ are the unperturbed 
eigenenergies of ffTOl) . The spontaneous decay rate is obtained applying Fermi's golden 
rule to the interaction Hamiltonian of each resonator with its harmonic bath. The 
transition rate from the coupled transmons plus resonators state |cr, Oi, O2) + |cr, Oi, 02)^^-' 
to the ground state \g, Oi, O2), is 



2 



2^^,!- ' (11) 

rv ^ 



where is the spontaneous emission rate of oscillator a and we considered single photon 
losses to each bath. The coupled transmons SWAP levels experience a Purcell induced 
spontaneous emission rate reduced with respect to the sum of the resonators spontaneous 
emission rates. For identical transmons in cavities with a lifetime I/kq, ~ 160 ns we 
estimate l/wa ~ 16/is (analogously to the transmon's relaxation time predicted in 
Ref. [13]), signaling a limitation to the optimized gate efficiency. We expect that the 
recently proposed Purcell filter or protection schemes [25] can be suitably extended to the 
considered two-qubit gate which is based on independent readout, possibly overcoming 
Purcell limitation. 



5. Conclusions 



In conclusion, we demonstrated optimization of a cQED entangling gate against any 
relevant 1// noise source while keeping the hardware simplicity of the fixed coupling 
and even if one qubit does not operate at optimal bias point. The estimated high 
performance of the gate signals the effective elimination of leading order effects of 1// 
noise. 

Our analysis included all the relevant noise sources acting during the entanglement 
generating operation in the considered architecture. We have shown that the proposed 
scheme is robust with respect to relaxation processes due to quantum noise and it is likely 
to foresee a design protected also from Purcell effect. Additional errors during readout 
may of course influence the overall gate fldelity of any speciflc implementation [9] . The 
responsible error sources need to be independently eliminated. However, the value of 
the optimal coupling is not affected by minimization of error sources acting before/after 
the coupled-qubits evolution. Similarly, for qubit-based quantum information |36j . 
optimization of single and two qubits quantum operations is a key requirement, even 
though the overall quantum processor will suffer from error sources in between quantum 
operations or at preparation/readout. 

Eliminating decoherence remains the biggest challenge for superconducting systems. 
Further optimization may require on one side suppression of higher order effects of 
1// noise, on the other limitation of relaxation due to quantum noise. Concerning 




Figure Al. Absolute value of the coherence {p^ (t)) / (0) in the presence of 1// 

flux noise on qubit 2 with Y,X2 = 10~^. The (thick) red line is the result of the SPA, 
the (thin) blue line is the numerical evaluation of the adiabatic approximation (|A.1|) 
for = ls~^, 7A/2 = 10^ s~^ as in [31], the dashed blue line is for 7^2 = Is^^, 
7A/2 = 10^ s"-'^: the smaller is the high-frequency cut-off 7m2, the closer is the 
SPA to the adiabatic approximation. Inset: zoom around the time range where 
/9+_(t)//9+_(0) sa e-^. Other parameters are Eca = 1 GHz, E'}^ = 30 GHz, 
Ecc ~ 1-68 X 10^ Hz {Ecc ~ 65 MHz). In the simulations we considered an ensemble 
of ~ ID'' random telegraph noise processes with switching rates distributed as cx I/7 in 
[7m2,7A/2]- The average is performed using 10^ realizations of the stochastic process. 



intrinsic noise sources, like those responsible for 1/f noise, material engineering at the 
microscopic scale may be required in the near future. "Passive" optimization startegies 
may be conveniently combined with "active" control tools, like dynamical decoupling 
protocols inspired to nuclear magnetic resonance which have been already applied to 
superconducting systems [311 [37]. On a longer time scale, imperfections in the coherent 
control might represent the ultimate limit to computer performance. 
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Appendix A. Validity regime of the Static Path Approximation 

The SPA is valid for times t < Thus it applies to the considered \/i — SWAP 

operation if = 7r/2cj+_ ^ 10^^ — 10~^GHz < l/'jMa- Since flux noise is the most 
relevant 1/f noise source in the considered setup, here we disregard critical current 
fluctuations and consider the recent noise figures reported in Ref. [21] • In that article 
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1// flux noise extends up to ~ 1 MHz, thus we can reasonably expect that the condition 
tE < '^/iMa is satisfied. 

Moreover, we numerically verified that the SPA is a valid approximation also for 
times t > 1/7m provided that 7m is smaller than the system oscillation frequency [29] . 
In figure RD we report the coherence between the states |±) in the SPA and the result 
of the numerical evaluation of the adiabatic approximation 

= p+_(0) j V[{x^{s)}] P[{x^{s)}] e-^Io^^-.-(M^)}) , (A.l) 

where uj+-{{xa{s)}) ^ Y.a + 1 Y.a,p fx'^dxp ^o,{s)xi3{s) and the derivatives are 

reported in Section [3] below Eq. (Q. In the figure we considered fiux noise on qubit 2, 
^2(5), distributed with = 10""^ and 7^2 = 1 s""*^, 7^/2 = 10® s~^ [31]. It is clearly seen 
that the SPA is a reasonable approximation up to times ~ 10^/7Af2- This legitimates 
the use of the SPA for the evaluation of the times reported in figure [31 The 

error with respect to an estimate based on the adiabatic approximation is of ~ 10% 
(inset of figure lAip . 
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